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Abstract 

A non-subtractive recipe of Casimir energy renormalization efficient in the presence of logarith- 
mically divergent terms is proposed. It is demonstrated that it can be applied even in such cases, 
when energy levels can be obtained only numerically whereas neither their asymptotical behavior, 
nor the analytical form of the corresponding spectral equation can be studied. The results of numer- 
ical calculations performed with this method are compared to those obtained by means of explicit 
' subtraction of divergent terms from energy. 
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1 Introduction 



Ever since Casimir [1] has obtained corrections to the energy of a macroscopic system due to vacuum 
fluctuations of quantized electromagnetic field in 1948 this effect has been intensively studied both from 
■ theoretical and experimental points of view. Nevertheless the calculation of Casimir energy except for the 

most simple problems involving free fields inside cavities with flat boundaries remains quite non-trivial 
yet. To realize this, one can recall a great number of papers devoted to a free field confined in the interior 
of a sphere [2]— [11]. Despite the fact that in this case spectral equations can be written out explicitly, the 
first analytical results for massive scalar field have been obtained only in [12]. The analogous problem 
for fermions has been solved in [13]. 

It should be stressed that the knowledge of analytical form of spectral equation has been crucial for 
the employment of the method proposed in [12] since it makes possible the transition from the sums 
containing the unknown energy levels to the integrals with the explicit integrands [14]. The main goal 
of this paper is to propose a method which can be applied to numerical calculation of Casimir energy 
in situations when these requirements are not met. Another goal is to demonstrate that it's possible 
to perform all the necessary calculations without employment of a rather standard trick [7], [15], which 
lets one overcome problems arising due to the presence of logarithmic divergency in Casimir energy of 
free massless fields inside spherical shells. Note that logarithmic divergency appears as a consequence of 
a curved surface bounding the shell and makes the energy renormalization ambiguous. The main idea 
of the mentioned trick is to consider the "inner" and "extra" problems together since their logarithmic 
divergencies cancel each other. 



The ambiguity of Casimir energy renormalization in the presence of logarithmic divergency is quite 
obvious. Indeed, in case of massless fields the energy of the system can be characterized by a single 
dimensional parameter L which is the linear size of the system. The regularization parameter a can 
be also chosen to have a dimension of length. In the absence of logarithmic divergency the "minimal 
subtraction" of singular terms is not only natural but also well grounded. Indeed, any divergent term in 
the expansion of regularized energy, which is proportional to oT s (s > 0) is inevitably proportional to 
L s , i.e. to the non-negative power of L. This makes it possible to normalize the final result at L = oo, 
where Casimir energy should become zero, and subtract all singular terms at the same time. After such 
subtraction the only remaining term in the limit a — > 0, which reads ca° /L, provides the final result. 

In the presence of logarithmic divergency the subtraction becomes ambiguous, since in order to renor- 
malize the term ca° log(a/L)/L, one should subtract ca°log(da/ L)/ L, where d is an arbitrary constant, 
which cannot be determined from the normalizing condition at L — > oo. 

2 Massive scalar field in ID 

To illustrate the main idea of technique under investigation let's study Casimir energy with the logarithmic 
divergency in the most trivial case, i.e. Casimir energy of massive scalar field on an interval of length L 
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with Dirichlct boundary conditions at the ends of the interval. It reads: 

1 1 

n— 1 n— 1 

First of all, the renormalization of Casimir energy is performed without intermediate subtraction of 
Minkowski vacuum contribution prescribed by "standard" approaches. The regularization of (1) requires 
the introduction of parameter a having a dimension of length, which stands in the argument of the cut-off 
function F(au) n ): 

oo 

£i:l = ^u n F(aw n ) (2) 

Z n=l 

Trivial considerations based on dimensional analysis lead to the following expression for the regularized 
Casimir energy 

L LP 1 

£gl ~ c_ 2 — +c^—+ c - + c x m 2 Llog(a/L) + ■■■ (3) 
or a 1 L 

It can be easily verified that for various cut-off functions such as F(x) = cxp(— x), F(x) — exp(— x 2 ), 
F(x) = cxp(— x 3 ), . . ., F{x) = exp(— a; 6 ), F{x) = cxp(— 2 cosh(x) + 2), . . ., identical c\ are obtained, while 
Co are different. The identity of c\ for different F(x) can be demonstrated by means of the following 
estimation for the sum giving rise to the logarithmic divergency: 
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F(aw n ) ~ — — log TV ~ — — \og(AxL/a) = 

Z UJ n 47T 47T 

■ \og(L/a) + — — log(Ax) , (4) 
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where N ~ AxL/a and Ax is a cut-off interval of F(x). 

Since any subtraction in the presence of logarithmic divergency is ambiguous this procedure should 
be excluded from consideration along with the logarithmic divergency itself. To achieve that one should 
calculate d\£^ s : 

d 2 A r a l^c ^+c x ^ + -.- (5) 

The obtained expression is regular in the limit a — > 0, so no subtraction is required. The knowledge of 
the function d^Scas, lets one reconstruct the required £cal unambiguously, since the initial conditions at 
L — > oo are well-known: both £cas and OlSccis should become zero in this limit. Note that while (5) 



doesn't describe the asymptotical behavior of dlScal at L — > oo, it demonstrates the disappearance of all 
singular terms in this quantity. Moreover the following integral approximation shows that d^Scal does 
really vanish for L — > oo: 

_^ oo _^ oo 

n—1 n— — oo 

k\ [ dx(L/Tr)Vx 2 + m 2 F(aVx 2 + m 2 ) - ^F(am) (6) 
2 J-oo 4 



3 Method of calculation in general case 

The proposed method can be generalized in such a way that it doesn't require the analytical expression 
for energy levels for its application. Suppose one has a set of energy levels of some spectrum uj n and 
the corresponding regularized Casimir energy contains the logarithmically divergent term. First of all it 
turns out to be possible to modify the initial expression for the Casimir energy by introduction of some 
parameter /i in such a way that 
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1 oo 

-J2V^T7F(a^4^) (7) 
1 „=i 

no longer contains the logarithmic divergency. For the massive scalar field on an interval /x obviously 
equals to the mass of the field. In less trivial three-dimensional cases with spherical symmetry /i is 
some parameter having a dimension of mass which characterizes the total coefficient by the logarithmic 
divergency with all values of angular momentum taken into account. 

The next step is to introduce an "additional mass" of the field M. and study the behavior of Casimir 
energy in the range from A4 = to A4 = oo. In fact it's convenient to deal with another parameter M, 
which is related to M. by means of M 2 = M 2 + ^ 2 , and study the "modified" Casimir energy 

oo 

S cas (M) = -J2 y/u*+M*-f (8) 

1 n=l 

as a function of M in the range from M = oo to M = fi. To carry this out, one should calculate 
numerically the following quantity in the specified range of M: 
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d 2 M (£ { c :l(M)/M) = 

I,, a - ,,2 V 

(9) 



n=l \ 



Note that in contrast to our previous considerations the dimensionless quantity 



has been substituted to the argument of F(x), so that a should be also taken dimensionless. 

An alternative interpretation of (9) follows from the observation that M acts just as an effective length 
L in the expression for Casimir energy. Indeed, (9) can be obtained from the initial expression as a result 
of the following sequence of transformations of the spectrum. At the first step u„ is transformed to 
ui' n = y/u) 2 — fx 2 / n, which is the dimensionless form of the spectrum with the subtracted effective mass. 
After that the scale transformation of the system x — > x{Mj \x) and oj' n — > = u' n /(M/fi) is performed. 
In the end the unit mass is "added" to the obtained spectrum: 



<-VK) 2 + i = ^ + i (io) 

The limit M — > oo obviously corresponds to the infinite size of the system, while for M = \i one obtains 
the initial spectrum divided by /U. 

It's easy to see that all divergent terms in (9) vanish. In the limit n — > oo one can make use of the 
following expansion 

HZ7i~^v^! + — +-■ (id 

M 2 M + 2^J^2 W 

The first term in this expansion gives rise to the sum which is free of logarithmic divergency owing to the 
definition of fi. Other divergencies are proportional to (M/a) 2 /M and (M/a) 1 /M and vanish when the 
second-order derivative is taken. The second term leads to the logarithmic divergency with the coefficient 
by it proportional to M 1 , which also vanishes while taking the second-order derivative. 

As a result one has the expression (9) regular in the limit a — > and the natural normalizing condition 
£ cas (M — > oo) = 0. The latter can be understood from two different points of view. On one hand the 
quantized field with the infinitely large mass should definitely have zero Casimir energy. On the other 
hand the Casimir energy in the limit of the infinite size of the system should become zero. Whichever 
interpretation is chosen, these two points let one reconstruct the required £ cas (M — /i) corresponding to 
the non-modified spectrum. 
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Note, that principally one could consider £ca S (M) instead of Seal {M) /M. However that would increase 
the order of derivative required to exclude all divergent terms what is undesirable from the point of view 
of real-time computations. 

The proposed method turns out to be efficient not only in the most trivial one-dimensional cases, but 
also in more realistic three-dimensional ones. However to employ it in three-dimensional case one should 
inevitably calculate the fourth-order derivative of the Casimir energy (9) since the main singular term, 
which is proportional to the volume of the system, reads c_ 4 L 3 /a 4 . It should be also noted that in this 
case the summations become more lengthy and sophisticated. Typically the final result turns out to be 
about 40 orders lower than the values of partial sums obtained in the process of its calculation. As a 
consequence, extra floating-point precision is required. 

4 Numerical results 

For scalar field on an interval [0; L] with L = 1 the spectrum reads 

Wn = ^ + « 2 (12) 
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Fig. 1: The Casimir energy of the massive scalar field on a unit interval as a function of the mass of the 
field. 

The result of straightforward application of the proposed technique with various cut-off functions 
such as F(x) — cxp(— x), F(x) — exp(— a; 2 ), F(x) = exp(— x 3 ), . . ., F(x) — exp(— x 6 ), F(x) = 
cxp(— 2cosh(x) + 2) is presented on Fig. 1. It has been shown that for each of these functions the 
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same result is obtained and what's more the precision of coincidence depends only on the number of 
energy levels taken into account and the number of right digits used in the realization of floating-point 
arithmetics as well. 

As to dependence of the Casimir energy on the mass of the field some important aspects should be 
stressed. First of all in the limit m-»0a well-known result for the massless scalar field is obtained. In 
the range of large values of m Casimir energy decreases exponentially as e ~ 2mL what could be expected 
from qualitative considerations. To summarize, the results obtained in this trivial case are in full agree- 
ment with the previously known results, obtained with the traditional technique, based on an explicit 
subtraction of divergent terms. 

To demonstrate how this approach can be employed in less trivial cases it makes sense to consider the 
massless scalar field inside of a spherical shell of radius R = 1. In this case the same set of cut-off functions 
has been exploited. The values of an effective mass fi — 0.1377 obtained with each of these functions 
coincide up to the first four digits. Consequently the precision of the obtained E cas — 3.790 • 10~ 3 
has the same order, what corresponds to about 200 s-levels taken into account while performing the 
calculations. In fact, the number of energy levels taken into account is directly affected by the range 
which the regularization parameter a used in the calculations belongs to. Therefore one can control 
precision of the final result simply changing the range of employed values of the regularization parameter. 

Note that in the framework of this approach not only Casimir energy of the massless scalar field 
(corresponding to M = 0) has been obtained but also Casimir energy for all possible mass values in the 
range from zero to the "effective" infinity. The dependence of the Casimir energy of the scalar field inside 
the sphere on the mass of the field is presented on Fig. 2. 




Fig. 2: The Casimir energy of the scalar field inside of a spherical shell of radius 1 obeying Dirichlet 
boundary conditions as a function of the mass of the field. 
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It should be stressed that on the contrary to the results obtained with the traditional subtractive 
technique in [12] our result doesn't contain logarithmical singularity at M. = what seems more rea- 
sonable from the physical point of view. The most likely explanation of this discrepancy is that the 
subtractive procedure may very well contain some arbitrariness, even if one puts the proper normalizing 
condition. As a result some function which has regular behavior at A4 — > oo (and thus doesn't violate 
the normalizing condition), but is singular at M — > may have been subtracted from the final result. 




Fig. 3: The Casimir energy of a massive scalar field in the whole space with Dirichlet boundary condition 
on a sphere of radius 1 as a function of the mass of the field. 

As has been pointed out in [12], [15] there is no argument at present which can remove this arbitrariness 
in case of a massless scalar field in the interior (or exterior) of a sphere. Therefore Casimir effect in the 
whole space with Dirichlet boundary conditions on the sphere is usually considered instead of interior 
alone (see the Introduction). It seems instructive to calculate the Casimir energy in this specific case 
employing our technique. 

In fact there arc two ways to proceed to take exterior into account. The first one deals with the 
continuous spectrum and requires that all the regularized sums be replaced with the appropriate integrals 
containing the energy levels density in their integrands. The second way lets one work with discrete 
spectrum all the time. To carry that out one should place the initial spherical shell into another sphere 
with the radius R ou t — kRi n where k > 1 and calculate the Casimir energy for the system bounded by 
the outer sphere with the boundary conditions on both of the spheres taken into account. For each finite 
k the spectrum remains discrete and the developed technique can be applied without any modification. 
The required result can be achieved in the limit k —* oo. In practice it turns out that for k > ko, where 
k is finite and depends on the required precision only, the result doesn't depend on k. It turns out that 
in the case under investigation the 4-digit precision of the final result can be achieved with ko — 10. 
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The final result of the calculations is presented on Fig. 3. Note that while the qualitative behavior 
of Casimir energy is the same as that obtained with methods employing explicit subtraction [15], there 
is no absolute coincidence. For example, for the massless field the result obtained with our method is 
£ C as(M = 0) = 0.0039 while direct subtraction leads to £ cas (M = 0) = 0.0028. Again, the remaining 
discrepancy can be related to the finite arbitrariness of the subtractive procedure. 

5 Conclusion 

To summarize, an efficient technique for numerical calculation of Casimir energy of quantized fields in 
the presence of logarithmical divergencies, owing to non-trivial boundary conditions, has been developed. 
The advantages of the proposed method are its ideological simplicity and universality which let one apply 
it to a wide range of problems in which numerical values for energy levels can be obtained with the 
sufficient precision. As has been demonstrated the results of its application to a number of problems 
appear to be quite reasonable, including the cases with curved boundaries. As to possible disadvantages 
they turn out to be mostly technical: the required floating-point precision to carry out all the necessary 
calculations in realistic cases turns out to be quite high (from three to four times higher than the one 
realized in the standard C double type). 
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